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ABSTRACT 

We present evidence that the spectroscopically identified bipolar jets of the pre-main sequence binary 
KH 15D (P=48.4d, e periastron separation --18 Ra, Ma=0.6Mo, Mb=0.7Mo) are a common product 

of the whole binary system, rather than being launched from either star individually. They may be launched 
from the innermost part of the circumbinary disk (CBD) or may result from the merging of two outflows driven 
by the individual stars. This evidence is based on high-resolution Ha and [OI]A6300 line profiles obtained 
during eclipse phases of this nearly edge-on system. The occultation of star A (the only currently visible star) 
by the disk strongly suppresses the stellar Ha and continuum emission and allows one to study the faint red- 
shifted and blueshifted emission components of the bipolar jets. The strongest evidence for jet production by 
the whole binary system comes from the observed radial velocity symmetry of the two jet components relative 
to the systemic velocity of the binary in combination with current accretion models from the CBD onto a binary 
system. 

Subject headings: binaries: spectroscopic - circumstellar matter - ISM: jets and outflows - stars: individual 
(KH 15D) 



1. INTRODUCTION 

The models proposed for driving jets and outflows from 
young stars can be roughly categorized as disk-driven out- 
flows (e.g. Pudritz et al. 2007), X-winds (e.g. Shang et al. 

2007) , and stellar wind-driven flows (e.g. Matt & Pudritz 

2008) . In all of these models the young star is surrounded by 
a circumstellar disk with an outer radius of about 10 - 100 AU 
and an inner radius of only a few stellar radii. So far it has not 
been investigated whether any of these models can be applied 
to a close binary system such as KH 15D, where the spectro- 
scopic signature of a bipolar outflow has been found in the 
[OI]A6300 line (Hamilton et al. 2003). A small-scale bipo- 
lar jet of --70 AU length in [SII]A6731 is also known for the 
close 19d TTS binary UZ TauE (Hirth et al. 1997) and the 
[OI]A6300 Une profile of the 15.8 d TTS binary DQ Tau of 
Huerta et al. (2005) suggest that this system also drives an 
outflow. In this Letter we discuss for the first time simple 
possible models for jet formation in close TTS binary systems 
and hope that our Letter will initiate detailed theoretical and 
observational studies of this challenging astrophysical prob- 
lem. 
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In recent years KH 15D (age'--2Myr, d=760pc) has been 
the subject of many detailed observational and theoretical 
studies (e.g. Winn et al. 2006; Herbst et al. 2008 and 
references therein). It first got attention for its unique 
photometric variability with deep and practically grey 
eclipses every 48.4 days that have become deeper and wider 
with each passing year (Hamilton et al. 2005; LeDuc et 
al. 2010). As discussed in these references the following 
model describes the available data. A nearly edge-on binary 
system is surrounded by a circumbinary disk (CBD) which is 
inclined to the binary plane by «10-20° (Winn et al. 2004; 
Chiang & Murray-Clay 2004). Since more than a decade 
the CBD has been totally occulting the orbit of star B and 
has been increasingly occulting that of star A. This is due 
to its precession on a ^lO-'yr time scale, which causes the 
occulting edge to cover the orbit of the binary. During the 
epoch of observations discussed here, only eclipses of star A 
are observed, which result from the disappearance of this star 
behind the occulting edge. Without this sharp edge acting 
as a "natural coronagraph" it would have been practically 
impossible to trace in Ha or [OI] the much fainter emission 
from the outflowing gas analyzed here. 
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2. OBSERVATIONS 

The high-resolution spectra analyzed here have been ob- 
tained during various observing runs with UVES at the VLT 
(R-50000, slit width 0.8-1.0") and with HIRES at Keck 
(R--70000, slit width 0.8"). A full description of the data and 
their reduction wiil be given in Hamilton et al. (2010). At the 
distance of KH 15D, 1 .0" corresponds to 760 AU. The observ- 
ing dates of the individual spectra are indicated in Figs. 1 and 
2. The continuum levels of all these spectra have been cali- 
brated onto a common flux scale by using the Cousins 1 band 
fluxes that were measured simultaneously (or nearly so) in our 
CCD imaging campaigns. One uncertainty in this procedure 
is the system's R-I color, since R band fluxes are requested for 
calibration of Ha and [OI]. The R-I color was calculated on 
the basis of the 2002/2003 data set (Hamilton et al. 2005). In 
the spectra presented here, any uncertainty will most probably 
affect only the two spectra taken at minimum light (i.e. from 
Dec. 14, 2001 and March 12, 2004), since only there is die 
R-I versus I correlation unreliable. 

3. JET PROPERTIES DERIVED FROM [OI]A 6300 AND Ha 

It is now generally accepted that the forbidden emission 
lines (FELs) of T Tauri stars (TTSs) are an important diagnos- 
tic for their outflows close to the star on scales of 20-40 AU 
(Hirfli et al. 1997; Hartigan et al. 1995; Eisloffel et al. 2000). 
The most important line for these studies is the [01] A 6300 
line, which is usually the strongest EEL in TTSs in the opti- 
cal and therefore most easily observable. For KH 15D, the 
rather small measured out of eclipse Wa([OI]A6300)=0.4- 
0.5 A (which would be ^0.2 A if both stars would be unoc- 
culted), compared to classical TTS (CTTS) jet sources, indi- 
cates a mass outflow rate that is ^10 times smaller than in 
other CTTSs (e.g. Hartigan et al. 1995). 

In principle the Ha line is a much better diagnostic for 
investigating outflows on very small scales (e.g. 0.1-lOAU 
from the source) because it is not quenched above a certain 
electron density like EELs. Note that the critical density of 
the [OI]A6300 line of w 10'' cm is among the highest of all 
FELs observable in the red spectral region of TTSs. How- 
ever, in most TTSs the very bright steUar (magnetospheric) 
Ha emission makes studies close to the star extremely diffi- 
cult. For KH 15D, however the natural coronagraphic mask 
(see Sect. 1) strongly suppresses the stellar Ha emission and 
allows to study the faint Ha emission from the outflow. 

In Fig. 1 we display all our nine UVESA^LT [OI]A6300 
line profiles. They were all taken with UVES since the HIRES 
spectra do not cover this line. In the top four "uneclipsed" or 
"partially eclipsed" profiles the superimposed stellar absorp- 
tion spectrum has been subtracted to minimize it's effects on 
the profile shape (see caption for details). The three 2001 pro- 
files have already been published by Hamilton et al. (2003). 
In Fig. 2 we show 5 Ha line profiles obtained with UVES 
and HIRES on various nights during eclipse, which are quite 
different compared to "uneclipsed" profiles (Hamilton et al. 
2003, 2010). In both figures the given velocity scale is rela- 
tive to the systemic velocity {whei = 18.6km s~^). 

The nine [OI] profiles displayed in Fig. 1 are very similar 
and only small variations are indicated. In consideration of 
the night sky and flux calibration uncertainties, and the addi- 
tional noise in the top four profiles due to the absorption spec- 
trum subtraction, it is unclear, whether the Vrad values and the 
relative fluxes (f) shown in Fig. 1 indicate any significant vari- 
ation. The only exception is the profile from Dec. 20, 2001, 



where the flux tropped by about a factor of 2 (see Sect. 5 for 
an interpretation). In most profiles in Fig. 1 two distinct (but 
partly blended) peaks at ~-22km s"' and ~+25km s"' are 
visible. On the basis of the three 2001 [OI] profiles Hamil- 
ton et al. (2003) interpreted them in terms of a bipolar jet in 
which the two jet components lie close to the plane of the sky 
resulting in relatively small radial velocities. The Ha line 
profiles obtained during eclipse strengthen this interpretation, 
since several of them show signatures of bipolar jets similar to 
the [OI] line profile. A particularly nice example of a double- 
peaked profile is the middle one in Fig. 2 from Dec. 16, 2003. 
When comparing the [01] line profiles with the Ha line pro- 
files it is evident that the former indicate lower velocities for 
the bipolar jets. This could be explained by the presence of 
a blended and unresolved low- velocity component, which is 
often observed in the FELs of CTTSs, but not in Ha . 

Compared to [01] the presented Ha line profiles indicate 
considerable variability, particularly in the red jet component 
(which is often much fainter than the blue one) and in the 
broad Ha line wings. The data from all of our observing runs 
(each run sampling only a small fraction of the orbit) indicate 
only smaU variations on a 1-3 d timescale for the blueshifted 
jet component and larger variations only over months to years. 
The variations in the broad Ha wings are best illustrated by 
the unusual profile at the top of Fig. 2, with two additional 
broad peaks at -^ztlOOkm s~' and wings extending up to 
±250 km s"' . Similarly broad but much fainter wings are vis- 
ible in all Ha profiles, particularly on the blueshifted side. 
The nature of these variable Ha wings will be discussed in 
detail by Hamilton et al. (2010) in terms of magnetospheric 
accretion using all available profiles. 

To estimate the fuU opening angle a of the jets and the an- 
gle, ijet, between the flow axis and the line of sight we only 
consider the Ha line profiles, because the [OI] profiles could 
be blended with an unresolved low-velocity component (see 
above). For this purpose we assume that the jets can be ap- 
proximated by a homogeneously emitting cone. We will use 
Vrad and FWHM values of the Ha profile from Dec. 16, 2003 
to derive ije, and a from eq. 5 of Mundt et al. (1990). In this 
profile the two peaks are at -30 km s"' (FWHM=39km s"') 
and +34 km s~' (FWHM=48km s~^). If one assumes a jet ve- 
locity, Vj„, of 200 (100) km s"' one derives ;y«=81°(71°) for 
|v™d|=32km s-i and a values of 13°(26°) and 16°(33°) for 
the redshifted and blueshifted jets, respectively. 

4. POSSIBLE JET LAUNCHING MODELS 

Before discussing possible models we like to remind the 
reader that [Ol]A6300 emission in TTSs with small-scale jets 
is usually extending over 10-30 AU (e.g. Mirth et al. 1997). 
We call this spatial extend L^m, which is the length over which 
most of the light from the jet is emitted. Since the KH 15D jet 
has probably a smaUer mass flux than typical CTTSs in 
[01] is more likely lOAU flian 30 AU (see also Sect. 5 for 
an independent estimate of Le,,,). For Vye,=100km s"' it takes 
180 d to reach 10 AU. If one adopts for KH 15D such values, 
even a highly modulated outflow with Pi„„=48.4d would not 
cause detectable [OI] flux variations, because the [OI] emis- 
sion would be averaged out over ~3.7 periods. Although Ha 
is probably formed at smaller length scales (^2-5.5 AU; see 
Sect. 5) it is unlikely that Lem(Ha )/y jei is much shorter than 
the binary period and therefore one would expect only a small 
period modulation of the Ha flux of the jet component. These 
considerations are fuUy supported by the small [OI] flux vari- 
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ations indicated by Fig. 1 and the lack of periodic [Ol] flux 
and profile variations in the CTTS binaries DQ Tau and UZ 
TauE (Basri et al. 1997, Huerta et al. 2005). 

Due to the high excentricity (~0.6) of the KH 15D bi- 
nary system the two stars approach as close as ~18Ra 
(Ra=1.3R0) during periastron, i.e. near periastron any in- 
dividual circumstellar disks must be much smaller than this 
value due to tidal truncation. We find it of great astrophys- 
ical interest that such a close binary system can launch jets, 
and it remains as a large theoretical challenge to model such 
a system, or even closer systems like UZ TauE. One im- 
portant guidance for any launching model of the KH 15D 
jets is the fact that the redshifted and blueshifted jet com- 
ponents in Ha and [OI] show a symmetric velocity shift 
relative to Vjj,,=+18.6km s~^ In the following we discuss 
two possible models which could explain this important re- 
sult. In the first model the jets are launched from the in- 
nermost part of the CBD. A problem with this idea is the 
large inferred inner radius of the CBD (^0.6 AU; Herbst et 
al. 2008), since at this radius the CBD is rotating with only 
Vrof ~45km s"^ (P=145d, 3:1 inner Lindblad resonance). 
Magneto-hydrodynamic simulations (e.g. Ouyed & Pudritz 
1997; Fendt 2006) for various magnetic field strengths and 
magnetic field configurations have shown that jet velocities 
of Vjet ~1-1.5 Vror are normal, although Vye;~2vror is not un- 
usual. This means that vy<.f=100km s~^ is not implausible and 
we adopt this value for the CBD launching model. Accord- 
ing to the above ijet calculations a vjet value of 100 km s~' 
would imply ijet=l\° which is roughly consistent with the 
incUnation angle of the binary system (jj„>,=83-92°) derived 
by Winn et al. (2006), in the context of a specific geomet- 
ric model of the occulting edge. Better consistency might be 
achieved by allowing for a possible tilt (by ~10°) between 
the binary plane and the innermost part of the CBD, provided 
this tilt is oriented roughly along the line of sight. According 
to this model binaries with shorter periods should on average 
have higher jet speeds. This is in accordance with the [OI] 
line profiles of UZ TauE and DQ Tau, in which the edge of 
their blue wing indicate Vjet >150km s~^ and >120km s"\ 
respectively (Huerta et al. 2005). 

Although the CBD jet launching idea is attractive we also 
consider the following second model, which produces higher 
jet speeds. This model requires that each of the two binary 
stars launches a magneto-hydrodynamically driven outflow 
with ~200km s~' from each of the circumstellar disks and 
that each of these two outflows is of similar mass flux and ve- 
locity. At some distance from the system these two outflows 
merge to form a common jet. Whether this rather speculative 
idea would work in the end requires detailed modeling. One 
complication is the complex and unstable interaction of the 
magnetospheres of the two outflows which should lead to re- 
connection of the magnetic field fines. Since in this model the 
line widths of the jets should be correlated with the radial or- 
bital motions of the individual circumstellar disks (with -Vrad 
amplitudes of ~±50km s~^), it has to be shown, whether this 
model would agree with the small measured velocity widths 
of the jet components (FWHM=40-50km s"'). This model 
predicts a periodic modulation of the outflow (see below), 
which would only be observable with an outflow tracer emit- 
ting very close to the source (i.e., Le„Jv ja <C Ptin)- 

According to the models of Gtinther and Kley (2002) bi- 
nary stars with similar masses should have similar mass ac- 
cretion rates Mncc from their CBDs. Due to the excentric orbit 



(e~0.6) Mace should be highly modulated with the binary pe- 
riod {Macc,maxl^acc,mm ^4-5) and highest near periastron pas- 
sage. Since in TTSs M outflow is highly correlated with Mace, 
one would expect for KH 15D similar Momjiow values for the 
two stars, which should in addition be highest near periastron. 
If, quite in contrast to theoretical expectations, only one of the 
two stars in the system would still drive the jets one would ex- 
pect a strongly modulated outflow with a higher fractions of 
the matter ejected near periastron. Since at this phase each of 
the binary members has a rather high Vrad ( ±50kms-' rela- 
tive to ^sys) it would be hard to explain why the two jet com- 
ponents are not symmetric about a different velocity, despite 
the fact that the recorded jet emission traces gas ejected over 
several orbital periods. These considerations make it rather 
unlikely than the jets of KH 15D are launched from only one 
of the two binary stars. 

5. VARIATIONS OF THE REDSHIFTED JET IN Ha 

The large variations in the "eclipsed" Ha line profiles and 
fluxes relative to the much more stable [Ol] line profiles in- 
dicates that these two lines cannot be formed in the same re- 
gion. Most probably the [OI] emission traces more distant 
and tenuous emission. This may mean that the [OI] emission 
is quenched in the inner and denser Ha emission region (i.e. 
Ng >>N^. „-,7 ~ 10'' cm"-') or that OI is largly ionized there. 
Whether any emission component in our spectra is modulated 
with the 48.4 d period cannot be answered with our data, be- 
cause of insufficient time sampling. A jet with 100 km s~^ 
would reach a distance of 2.7 AU within 48.4 d and, as already 
stated above, any periodic modulation in Ha would only be 
observable if ]^em in Ha is much shorter than this value. The 
shortest time scale variations of the red jet component have 
been observed between Dec. 14 and 20, 2001. The large flux 
increase of this component by a factor of 4 within 6 d cannot 
be explained by the flux correction uncertainties mentioned 
in Sect. 2. Perhaps they are related to the unique "accretion 
event" which happened on Dec. 20, 2001 (see Hamilton et al. 
2010). This could have produced sufficient UV radiation to 
excite more distant neutral sections of the jet causing a corre- 
sponding increase in Ha on a small time scale. The increased 
UV flux might have also ionized OI causing the smaller [OI] 
line flux observed on Dec. 20, 2001. 

We believe the large variations seen primarly in the red- 
shifted Ha jet component are mainly a result of obscuration 
by the disk combined with variation of l^em- Due to varia- 
tions in the outflow rate, or the above mentioned variations 
in UV excitation, L^.,,, for both the blueshifted and redshifted 
jet will change with time. If the redshifted jet is obscured 
by the circumstellar disk over a considerable fraction, \^obsc'> 
of its emitting section, then any changes in L(,„ will cause 
corresponding variations. If we assume that the outer disk 
radius, Rdisk, is at 5 AU (see Chiang & Murray-Clay 2004) 
and that the disk plane is inclined towards the line of sight 
by a=15° then Lo/„c=sin(a) Rdisk '^l.SAU. If one assumes 
Lem(blue)=Le„,(red) and the jet emission rate per unit length 
is constant, one can derive from the observed flux ratio of 
f(blue)/f(red)=f;/fr the ratio Le„,/Lo/„c=(fi,/fr)/((ffc/fr)-l), be- 
cause fi/fr=Lgm(blue)/(Le„,(red)-Lo|i,je). With the measured 
fhlfr values of ^^1.3-3.2 (and ignoring the unusual profile of 
Dec. 20, 2001) one derives for Ha Um ~2-5.5 AU. For [OI] 
the average ratio f^/f^ is 1.14. With the above equation one 
derives L(,m([OI])~10AU, in good accordance with the val- 
ues derived for other TTSs by spectro-astrometric methods. 
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Vrad ~ ^system [km/s] 

Figure 1. VLT [OI]A 6300 line profiles: the top four colored ones are taken out of eclipse (or at partial eclipse), while the five bottom ones are taken during 
echpse. In order to minimize profile changes due to the stellar absoiption spectrum, we have subtracted from the top four profiles a properly scaled and broadened 
spectrum of the K5 dwarf HD 36003. The number indicated on the right hand side (below each profile) gives the distance between star A and the occulting edge 
in radii of star A (for details see Herbst et al. 2008). The zero flux levels of the profiles are indicated on the left hand side (dashed colored lines for the four top 
profiles, black lines for the five bottom ones). The velocities of the jet components are indicated in km s"' and are based on Gaussian component fitting. The 
given f-values indicate the relative fluxes of the jet components, with the top right-hand component set to f=l. 
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Vrad ~ ^system [km/s] 

Figure 2. Sample of Ha line profiles obtained during eclipse at various observing runs at the VLT and KECK. The VLT [OI]A 63001ine profile from Dec. 18, 
2004 is shown for comparison (dotted red line). For details see the caption of Fig. 1 . 



